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Abstract. 

While the rising flux tube paradigm is an elegant theory, its basic assumptions, thin flux 
tubes at the bottom of the convection zone with field strengths two orders of magnitude above 
equipartition, remain numerically unverified at best. As such, in recent years the idea of a forma- 
tion of sunspots near the top of the convection zone has generated some interest. The presence 
of turbulence can strongly enhance diffusive transport mechanisms, leading to an effective trans- 
port coefficient formalism in the mean-field formulation. The question is what happens to these 
coefficients when the turbulence becomes anisotropic due to a strong large-scale mean magnetic 
field. It has been noted in the past that this anisotropy can also lead to highly non-diffusive 
behaviour. In the present work we investigate the formation of large-scale magnetic structures as 
a result of a negative contribution of turbulence to the large-scale effective magnetic pressure in 
the presence of stratification. In direct numerical simulations of forced turbulence in a stratified 
box, we verify the existence of this effect. This phenomenon can cause formation of large-scale 
magnetic structures even from initially uniform large-scale magnetic field. 
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1. Introduction 

The standard explanation for the appearance of strong magnetic fields at the solar 
surface involves the coherent rise of a tachocline-generated magnetic flux tube through 
the solar convection zone. Flux tube emergence simulations do give very promising results 
(jRempel et al.ll2009l) , the paradigm looks elegant and is very textbook friendly, but some 
of its assumptions are problematic: the integrity of flux tubes, their rise, and even their 
very existence. 

So far, numerical sim ulations have failed to produce the assumed thin magnetic flux 
tubes in the tachoclinc (Catta neo et al.1 2006; P arker 2009). Magnetic buoyancy as the 
driving transpo rt mechanism through the convection zo ne can be dominated by down- 
ward pumping (jNordlund et al.lll992l ; lTobias et al.lll998l ). For the tubes to r emain inta ct, 
strongly super-equipartition field strengths and strong twists are required ( Fanll2001 ). 

Thus, it seems not unreasonable to explore alternative mechanisms of formation of 
strong magne tic fields at the sol ar surface. A number of positive arguments has been 
put forward ( Brandenburgl 2005 ) in favour of formation of sunspots from local flux 
concentrations near the surface. In mean-field models, such magnetic instabilities have 
been produced by introduci ng a magnetic dependence of the thermal eddy diffusiv ity 
(Kitchatinov & Mazur l2000h and the viscous stress tensor (jBrandenburg et al]|2010ah . 

Turbulence generally is associated with enhanced transport effects. However, it also 
exhibits non-diffusive behaviour, generating magnetic field on much larger scales than 
the driving scale. An example is turbulent dynamos that are able to produce large-scale 
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magnetic fields ( Brandenburg fc Subramanianll2005 ). In a stratified layer, magnetic fields 



tend to become buoyantly unstable: assuming a constant temperature across the magnetic 
flux tube, pressure balance implies lower densities in regions of stronger magnetic field. 
Now one can wonder what the role of turbulent pressure is in such an equilibrium if 
the magnetic structures extend over several turbulent eddies. The turbulent pressure 
associated with the convective fluid motions is certainly not negligible and is strongly 
affected by the background magnetic field. The latter can be seen by evaluating the total 
turbulent dynamic pressure, here given for the isotropic case: 

ftwb = \pu 2 + \b 2 /no, 

here it and b are the velocity and magnetic fluctuations, respectively, /iq the vacuum 
permeability and p the fluid d ensity. Overbars indicate e nsemble averaging. As shown in 
direct numerical simulations ([Brandenburg et al. 2010a, hereafter referred to as BKR), 



in forced turbulence with imposed uniform large-scale magnetic field, the total turbulent 
energy is approximately conserved in this parameter regime. 

\pu 2 + ±b 2 /(i = E tot w const. 

As a result, one finds a reversed feedback from the magnetic fluctuations on the turbulent 
pressure: 

-Pturb = -gftV/XQ + 2Etot/3 

( Kleeorin et al. 19901 : Rogachevskii fc Kleeorinll2007 , hereafter referred to as RK07). One 



can see that the effective mean magnetic pressure force is reduced and can, in a certain 
parameter range, be reversed. This effect would then counteract the aforementioned 
buoyancy instability in the presence of turbulence. RK07 have suggested that the reversed 
feedback instability could lead to the formation of magnetic flux concentrations near the 
solar surface. 

Mean-field magnetohydromagnetic simulations by BKR confirmed the basic pheno- 
menon of magnetic flux concentration by the effect of turbulence on the mean Lorentz 
force and for sufficient stratification, a linear instability was found. Direct numerical 
simulations (DNS) by BKR confirmed the reversed feedback phenomenology, but did not 
address the e ffect of stratification . This i s one of the important additions of the more 
recent work of iBrandenburg et al. ( 2010b ). of which we report here the main highlights. 



2. DNS model and analysis 

DNS of forced turbulence were performed in a cubic computational domain of size L 3 . 
For an isothermal equation of state and a constant vertical gravitational acceleration g, 
one finds an exponentially stratified density: 

p = p Q exp(-z/H p ) , 

where H p — c 2 s /g is the constant density scale height, c s is the isothermal sound speed 
and po a normalisation factor. We choose k\H p = 1, where k\ — 2tt/L, the smallest 
wavenumber. The density contrast then corresponds to exp27r s» 535. 

We solve the equations of compressible magneto-hydrodynamics in the form 

P^r = J x B - clV Inp + V • (2^S) + p(f + g), 



— = U x B + f]V 2 A, 
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where v and r\ are kinematic viscosity and magnetic diffusivity, respectively. Furthermore, 
B = B + V x A is the magnetic field consisting of a uniform mean field, Bo — (0, B , 0), 
and a nonuniform part that is represented in terms of the magnetic vector potential A, 
J = V x B/fi is the current density, and Sy = \(Ui,j + — \hj V ■ is the traceless 
rate of strain tensor, where commas denote partial differentiation. The turbulence is 
driven with a forcing function / that consists of random plane non-polarized waves 
with an average wavenumber kf — 5 k\ . The forcing strength is arranged such that the 

1/2 

turbulent rms velocity, it rms = (u 2 ) , is around 0.1 c s . This value is small enough so 
that compressibility effects are weak. 

In order to characterize our simulations, a set of dimensionless parameters is defined. 
The Reynolds number is given by Re = u rms /iskf and is of the order 120 in the simula- 
tions. The magnetic Prandtl number is Ptm = v/rj, although in reality this parameter 
is much smaller than unity; we choose here Ptm = 0.5 — 8 in order to achieve higher 
values for the magnetic Reynolds number. The equipartition field strength B cq is defined 
as a function of z and the imposed fields are normalised against -B eq o, the equipartition 
strength in the middle of the domain: 

1/2 



B cq 



(z) = (a*0P" 2 ) , B cq0 = (^ Po) 1/2 "r 



The boundary conditions are stress-free perfect conductors at the top and bottom 
of the domain, periodicity in the horizontal direction. The simulations are performed 
with the Pencil Code0, which uses sixth-o rder explicit finite differences in space and a 
third-order accurate time stepping method ([Brandenburg fc Doblerll2002i ). 

The contribution to the mean momentum density by the fluctuations is 

Ilfj = py~u] + ^'iJ' 1 - h,b r 

where the overbars now indicate horizontal averages. The influence of the mean magnetic 
field can be found by subtracting the contributions that are present in the absence of a 
uniform background magnetic field and can be modelled by the following ansatz (RK07) 

n fj B ~ n f/ = ~ {\ 5 ij1p + e i e 3<le^ B 2 + q s BiBj, 

which will then appear in the effective mean Lorentz force 

pFf = -v, (V,77 ;! • 77,77, ■ tt;; /; nf/ 

Thus the total effective magnetic pressure of the mean field is given by 1 — q p (i?) B" As 
B fa (0,75,0), and assuming no small-scale dynamo, we can determine q p from 



( u l ~ U 0x) 



\V-bl = -\q p B\ 



3. Results 

We compare the effective magnetic pressure of the mean field with the turbulent kinetic 
energy density (Fig. [T^,) and see that the resulting contribution of stratified turbulence 

f http : / /pencil-code . googlecode . com 
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Figure 1. Normalized effective mean magnetic pressure for Bo — 0.fB cq o, Bo — 0.2B cq o, and 
Bo = 0.5.BcqO using Re = 120. Left: as a function of depth Right: as function of the local value 
of the ratio of Bo/ B cq (z). Note that the curves for the different imposed field strengths collapse 
onto a single dependence a nd agree very well with the fit by BKR (dashed line). Adapted from 
iBrandenburg et aU (|2010bl V 
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Figure 2. Left: Normalized effective mean magnetic pressure as a function of Bo / B cq (z) for 
Bo = 0.1-BcqO, and varying magnetic Prandtl number, Pr a/ = 0.05 (black), 2 (blue), 4 (red) 
and 8 (yellow). Right: Visualisation of B y — Bo on the periphery of the computational domain 
Bo = 0.1B cq o, and Pr M = 2. 

to the effective magnetic pressure is negative in a large part of the domain. Plotting the 
ratio of effective magnetic pressure to kinetic energy density as a function of the imposed 
uniform horizontal magnetic field divided by the equipartition field strength (Fig. [TJ)) , 
collapses the observations from simulations with different imposed field strengths into 
a single dependence of q p on B y /B eq . This result agrees with analytical calculations by 
RK07 and a fit based on simulations by BKR. The effect is fairly robust under an increase 
of the magnetic Prandtl number (Fig. [5^). While the reduced effective magnetic pressure 
is observed, the formation of local magnetic field concentrations as observed in mean-field 
simulations, has not yet been found in DNS (Fig. [2b). 

4. Discussion 

The DNS have shown that for an isothermal atmosphere with strong density strati- 
fication, the turbulent pressure is decreased due to a negative feedback from magnetic 
fluctuation generation, resulting in a negative effective mean magnetic pressure. The de- 
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pendence on the ratio of imposed field to local equipartition field agrees with results 
obtained from analytic theory (RK07) and direct numerical simulations (BKR). The 
results are robust when changing the strength of the imposed field and the magnetic 
Prandtl number. 

However, the simulations do not show any obvious signs of a large-scale instability that 
should result in magnetic flux concentrations, as was expected from mean-field calcula- 
tions. A possible explanation for this discrepancy could be the simplic i ty of the ansatz 



for the effective Lorentz force. Indeed, early simulations of iTao et al.l (|1998l ) produced 
clear signs of flux separation into magnetised and unmagnetised regions in convection 
simulations at large aspect ratio. An increase of scale separation might alleviate the effect 
of higher order terms. Future work should incorporate this increase as well as a wider 
scan of the parameter regime and ultimately the inclusion of radiative transfer, as was 



done in simulations of Kitiashvili et al. ( 201Clt ) , which showed flux concentrations in the 



presence of a vertical field. 

Acknowledgements 

We acknowledge the allocation of computing resources provided by the Swedish Na- 
tional Allocations Committee at the Center for Parallel Computers at the Royal Institute 
of Technology in Stockholm and the National Supercomputer Centers in Linkoping as 
well as the Norwegian National Allocations Committee at the Bergen Center for Com- 
putational Science. This work was supported in part by the European Research Council 
under the AstroDyn Research Project No. 227952 and the Swedish Research Council 
Grant No. 621-2007-4064. NK and IR thank NORDITA for hospitality and support dur- 
ing their visits. 



References 

Brandenburg, A. 2005, ApJ, 625, 539 

Brandenburg, A., & Dobler, W. 2002, Comp. Phys. Comm., 147, 471 
Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 

Brandenburg, A., Kleeorin, N., & Rogachevskii, I. 2010, Astron. Nadir., 331, 5 (BKR) 
Brandenburg, A., Kernel, K., Kleeorin, N., & Rogachevskii, I. 2010. larXiv:1005.5700l 
Cattaneo, F., Brummell, N. H., Cline, K. S. 2006, MNRAS, 365, 727 
Fan, Y. 2001, ApJ, 546, 509 

Kitchatinov, L.L., & Mazur, M.V. 2000, Solar Phys., 191, 325 

Kitiashvili, I. N., Kosovichev, A. G., Wray, A. A., & Mansour, N. N. 2010, ApJ, 719, 307 
Kleeorin, N.I., Rogachevskii, I.V., & Ruzmaikin, A. A. 1990, Sov. Phys. JETP, 70, 878 
Nordlund, A., Brandenburg, A., Jennings, R. L., Rieutord, M., Ruokolainen, J., Stein, R. F., & 

Tuominen, I. 1992, ApJ, 392, 647 
Parker, E. N. 2009, Spa. Sci. Rev., 144, 15 
Rempel, M., Schiissler, M., & Knolker, M. 2009, ApJ, 691, 640 
Rogachevskii, I., & Kleeorin, N. 2007, Phys. Rev. E, 76, 056307(RK07) 
Tao, L., Weiss, N. O., Brownjohn, D. P., & Proctor, M. R. E. 1998, ApJ, 496, L39 
Tobias, S. M., Brummell, N. H., Chine, T. L., & Toomre, J. 1998, ApJ, 502, L177 



